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Abstract 

We show that, during false vacuum inflation, a primordial magnetic field can be 
created, sufficiently strong to seed the galactic dynamo and generate the observed 
galactic magnetic fields. Considering the inflaton dominated regime, our field is pro- 
duced by the Higgs-field gradients, resulting from a grand unified phase transition. 
The evolution of the field is followed from its creation through to the epoch of struc- 
ture formation, subject to the relevant constraints. We find that it is possible to 
create a magnetic field of sufficient magnitude, provided the phase transition occurs 
during the final 5 e-foldings of the inflationary period. 
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1 Introduction 



One of the most exciting astrophysical consequences of phase transitions in the early uni- 
verse is the possible creation of primordial magnetic fields. 

The existence of a primordial magnetic field could have significant effect on various 
astrophysical processes. Indeed, large scale magnetic fields are important in intercluster 
gas or rich clusters of galaxies, in QSO's and active galactic nuclei. The existence of a pri- 
mordial field could influence the galaxy formation process and would play a very important 
role in the resulting galactic spins 0. A primordial field would also have an important 
effect on the fragmentation process of large scale structure and of the protogalaxies (by 
modifying the Jeans mass) and on the formation of Population III stars |J . But the most 
important consequence of the existence of a primordial magnetic field is that it can seed 
the observed galactic magnetic field. 

The galactic field is also very important to the astrophysics of the galaxy. It influences 
the dynamics of the galaxy, the star formation process (by transferring angular momentum 
away from protostellar clouds jj, Q and by affecting the initial mass function of the star 
formation process ||), the dynamics of compact stars (white dwarfs, neutron stars and 
black holes) and the confining of cosmic rays, to name but some. 

It is widely accepted that the galactic magnetic fields are generated through a dynamo 
mechanism usually referred to as the galactic dynamo, for which, though, there is no 
consistent mathematical model yet The basic idea of the dynamo mechanism is that 
a weak seed field could be amplified by the turbulent motion of ionised gas, which follows 
the differential rotation of the galaxy ||. The growth of the field is exponential and, 
thus, its strength can be increased several orders of magnitude in only a few e-foldings of 
amplification. 

The currently observed magnetic field of the Milky Way and of nearby galaxies is of 
the order of a fiGauss. If the e-folding time is no more than the galactic rotation period, 
~ 10 8 yrs then, considering the galactic age, ~ 10 10 yrs, the seed field needed to produce 
a field of the observed value is about ~ 10~ 19 Gauss ||, on a comoving scale of a 
protogalaxy. 

Although, it has been argued by many authors that the seed field could be produced 
from stars via stellar winds or supernovae and other explosions 0, there is evidence that 
suggests that the seed field is more likely to be truly primordial. For example, the observed 
field of the Milky Way does not change sign with z (z being the galactic altitude) as it 
would if it was produced by the stars of the galactic disk j3|. 

Various attempts have been made to produce a primordial field in the early universe. 
A thorough investigation of the issue was attempted by Turner and Widrow JF), who 
incorporated inflation and created the field by explicitly breaking the conformal invariance 
of electromagnetism. This was done in a number of ways, such as coupling the photon to 
gravity through RA 2 and RF 2 terms (R being the curvature, A being the photon field and 
F being the electromagnetic field strength), or with a scalar field 0, like the axion, through 
a term of the form (f)F 2 . It was, thus, shown that satisfactory results could be obtained only 
at the expense of gauge invariance. Garretson et. al. || have generalized the effort of [0] by 
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coupling the photon to an arbitrary pseudo-Goldstone boson, rather than the QCD axion. 
They have showed however that, in all cases considered, it was impossible to generate a 
primordial magnetic field of any astrophysical importance. Breaking the electromagnetic 
conformal invariance during inflation was a mechanism used also by a number of other 
authors, such as Ratra |9| and Dolgov |IJ. Ratra has been successful in generating an 



adequately intense magnetic field. The field was generated by coupling the field strength 
with a scalar field $ (the dilaton) through a term of the form e*F 2 . Dolgov, however, did 
not introduce any extra coupling but considered photon production by external gravity by 
the quantum conformal anomaly. He produced a field of enough strength, but only in the 
case of a large number (over 30) of light charged bosons. 

Another, more successful direction was using a phase transition for the creation of a 
primordial field. An early effort was made by Hogan |2|], who considered the possibility 
of turbulence arising during the QCD transition. His treatment, though, was based on a 
number of assumptions concerning equipartition of energy, which are of questionable valid- 
ity. Much later Vachaspati 0] proposed a mechanism to produce a marginally sufficient 
magnetic field during the electroweak transition. This has also been addressed, a bit more 
successfully, by Enqvist and Olesen [ P2]| . The later have also considered a phase transition 
to a new, ferromagnetic ground state of the vacuum, which could also produce an ade- 
quately strong magnetic field [13| . Finally, the literature contains a number of other, more 
exotic mechanisms (such as, for example, the creation of a primordial magnetic field by 



the turbulent motion of infalling matter into wakes in the wiggly string scenario [[14]]). In 



most of the cases, though, the achieved field appeared to be too weak to seed the galactic 
dynamo. 

In this paper we examine the production of a primordial magnetic field during false 
vacuum inflation. In false vacuum inflation a phase transition can occur during the infla- 
tionary period. As shown by Vachaspati [|11|] , the existence of a horizon could result in the 
production of a magnetic field at the phase transition. Although our model incorporates 
the benefits of inflation, it does not require the breaking of the gauge or even the conformal 
invariance of electromagnetism. Neither does it involve the addition of any extra couplings 
between fields through the inclusion of peculiar, "by hand" terms in the lagrangian. Our 
magnetic field is produced by the dynamic features of the GUT Higgs-field, independently 
of specific GUT models. No additional fields are introduced in the problem and the results 
cover the most general case. By considering a GUT phase transition we find that, for 
some parameter space, the magnetic field produced is of enough strength to seed a galactic 
dynamo mechanism at the epoch of structure and galaxy formation. 

Our results take into account any constraints imposed on the field during its evolution 
until the epoch of structure and galaxy formation. Before examining the behaviour of the 
magnetic field, we give a detailed description of the model and of the mechanism, through 
which the original primordial magnetic field is created. 
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2 Creation of the primordial magnetic field 



Vachaspati |TTJ suggested that the existence of a horizon would result in the creation of a 



primordial magnetic field at a phase transition in the early universe. 

Consider a non abelian group G. The field strength of the gauge fields is, 

Hpv = 9fj,X^ - d u X^ - QoffcX^Xl (1) 

where fg c are the structure constants of G and g is the gauge coupling. If the symmetry 
of the gauge group G is broken, leaving a residual symmetry corresponding to a subgroup 
H of G, then the gauge fields of the residual symmetry are given by, 

yb = u ab x a ^ 

where u ab is a unitary matrix specifying the directions of the generators of the unbroken 
symmetry. Vachaspati argued that if the symmetry was spontaneously broken then the 
VEV of the Higgs-field ip a would have been uncorrelated on superhorizon scales^], and 
hence, could not be "aligned" throughout all space with a gauge transformation. Therefore, 
the gradients of the Higgs-field would, in general, be non zero. Due to the coupling through 
the covariant derivative, 

D^) = (d,-igrXm (3) 

where r a are the generators of the residual symmetry and g is the gauge coupling, the 
corresponding field strength G is non-zero. This can be seen explicitly by using the 



gauge invariant generalization of 'tHooft Jl6] 



G\ v = u ah [H; v -g-^- 2 r cd D^ c D^ d ] 

= 9/^^y;-^yv/ c fl rf w d (4) 

where \i is the scale of the symmetry breaking. So even if the gauge field Y b can be gauged 
away the field strength is still non-zero: 

G\ v = -— 2 u ab f a cd d,rd^ d (5) 

Vachaspati applied the above in the case of the electroweak phase transition, taking 
G^ v to be the field strength of electromagnetism. 

In this paper we will follow a similar reasoning but for a GUT phase transition. G is 
now the GUT symmetry group and the residual symmetry is the electroweak. In order to 
get to electromagnetism we need to consider also the final electroweak phase transition. In 
analogy with the above the electromagnetic gauge potential is given by, 

A„ = v b Y b = sin 9 w n a W^ + cos 9 W B* (6) 



1 More precisely, on scales larger than the correlation length. 



3 



where v b is a unit vector specifying the direction of the unbroken symmetry U(l) em gen- 
erator, n a are the SU{2) generators of the electroweak group SU(2) x U(1)y, are the 
SU{2) gauge fields, is the U(l)y gauge field and 9w is the Weinberg angle. From (|) 
it easy to see that, 

v b = (sin 9 w n a , cos W ) (7) 
Y* = (W°BZ) (8) 



with b = 1, . . . , 4 and a = 1, . . . , 3 



The contribution to the electromagnetic field stength F^ from the GUT transition is, 
therefore, 



where v b G b is, 



F,u = v b G b (9) 



v b G b = sin 6 w n a G" + cos W G Y (10) 



B,\ ~ \F%,\ ~ |GW ~ -^(^(V)) 2 (12) 



with 6 = 1, . . . , 4, a = 1, . . . , 3 and = G^. 

The magnetic field produced by the GUT phase transition is, therefore, 

B, = \e^xF vX (11) 

Thus, for an order of magnitude estimate, (||), (H) and ( |I"l~D suggest, 

1 

since v a , u ab and /" d are of unit magnitude. As far as the Higgs-field gradients are con- 
cerned, on dimensional grounds, we have, 

*>M - 1 (13) 

where £ is the correlation length of the Higgs-field configuration. 

Also, at the GUT scale: 4:irg~ 2 ~ 40 =>- g^ 1 ~ 1 . Therefore, a dimensional estimate for 
the magnetic field is, 

B = iB„i~r 2 (14) 



3 False Vacuum Inflation 

In this section we review false vacuum inflation, a popular model of inflation corresponding 
to extensive literature (Hl7|...p3"|). 

In this model the inflaton field <fi rolls down its potential towards the minimum, which 
does not correspond to the true vacuum, but is instead a false vacuum state. There are 
two distinct and quite different kinds of false vacuum inflation, depending on whether the 
energy density is dominated by the false vacuum energy density or by the potential energy 



density of the inflaton field. Full details are given by Copeland et al in . Unlike them, 
we concentrate on the inflaton dominated case, as in |18"l|...|p2||. In this case, the phase 
transition does not lead to the end of inflation as it does in the vacuum dominated case. 
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3.1 The model 



In this model the energy density is dominated by a potential with two scalar fields; the 
inflaton field (f) and the Higgs-field ip. The later is responsible for the phase transition. We 
should emphasise here that the Higgs-field considered does not correspond to a specific 
GUT model and can have several components without this affecting the following analysis 

We take the form of the potential to be, 

V(<f>, = \\^ 2 - + + (15) 

The phase transition takes place at <fi = 0o> where, 

0o = ^ 2 (16) 
This gives the effective scale of the symmetry breaking, 

/4 = ^ - 1) (17) 

Without loss of generality we assume that <fi is initially positive and rolls down the 
potential in such a way that <fi < 0. If there is sufficient inflation before the phase transition 
and A ^> A' the Higgs-field will have rolled to the minimum of its potential ip = before 
the inflaton falls to its critical value cj) . So, when (f> > </> , 

n0,O) = iA// + imV (18) 
In the slow roll approximation the dynamics of inflation are governed by the equations, 



H 2 - {^V (19) 

3H4> ~ -V (20) 

where the prime and the dot denote derivatives with respect to <fi and time respectively, 
H = a/a is the Hubble parameter, a is the scale factor of the Universe, and m v \ is the 
Planck mass (m p i = 1.22 x 10 19 GeV). 
Thus, 

From (21) the number of e-foldings of expansion, which occur between the values <p\ 
and 2 of the inflaton field, is given by, 

ao 8% V 

N faM=**-r*-ZTj. ( 22 ) 
a\ m pl j (fix v 
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3.2 The inflaton dominated regime 

This is the case we are going to be interested in, since the back reaction of the Higgs-field 
ip to the inflaton field <fi is negligible and so the phase transition does not cause the end 
of inflation |p3| . If the opposite is true and inflation ends at the phase transition, then 
the effects of the transition are not too different from the usual, thermal phase transitions 
studied in the literature. 

In the inflaton dominated case the energy density of the inflaton field in (|i~8|) is much 
larger than the false vacuum energy density. Therefore, 

V{cf>) ~ im 2 2 (23) 



which is identical with chaotic inflation. Inflation ends at ~ <p £ , where, 

i _ m p i 



'Air 

To ensure inflaton domination until the end of inflation then, 



(24) 



i 



» W =► « 1 (25) 

2 £ 4 m 2 pl m 2 

This is the condition for inflaton domination. At this point it should be mentioned that 
if the inflaton domination condition is strongly valid the dynamics of inflation are not 
seriously affected by the phase transition, provided that ip falls rapidly to its VEV. Thus, 
the equations (|H)D, Q2DD and ([H]) can be used throughout the duration of inflation. 

Thus, for the number of e-foldings between the phase transition and the end of inflation, 
in the inflaton dominated case, we obtain, 

N = N(4> ,0) = (26) 
m P i 

where we have used <fi 3> 4> e . 

Finally, for the roll down of the inflaton field, using (|19|) , (f2~0|) and (|23|) we obtain, 



(27) 



4 The correlation length 

4.1 Evolution of the correlation length 

Through the use of the uncertainty principle, we can estimate the range of any interaction. 
Therefore, the physical correlation length for the Higgs-field ip is, 

t = (28) 
\m H \ 
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where m# = m#(t) is the mass of the Higgs particle, 



m # = — q^2 — = 3 ^ + V (29) 

Before the phase transition ip = and thus, m# = A|//^|. Therefore, the physical corre- 
lation length is, 

However, this is not valid as we approach the phase transition, <p — > O - As £ = 1 the 
correlation length grows linearly with time as shown in Figure 1. 
Define, 

T = t-t (31) 
where to is the time the transition occurs. 



If at r = t h , 

d , 1 



dr \vriH 

itic 

£o at the time of the transition is 



H 



1 (32) 



Then, from th until the transition r = 0, = 1 and, therefore, the correlation length 



£o = fn - th (33) 
where is the correlation length at r^-, and in the linear regime, 

e(r) = eo + r (34) 

The linear growth of £ continues until it hits the declining slope of (Figure 1). 

From then on the correlation length is given again by (|28|). 

The phase transition we are considering is not triggered by temperature fall, but by 
the roll down of the inflaton field. In that sense it is not a thermal phase transition. Also, 
since it occurs during inflation, the Universe is in a supercooled state with temperature 
T « and so there are no Ginzburg phenomena. However, the configuration of the Higgs- 
field ip does not freeze at the moment of the phase transition due to long wave quantum 
fluctuations that dominate the Higgs-field evolution immediately after the transition. 

The long wave fluctuations of ip are determined by the behaviour of the Higgs-field 
mass (P9|). Immediately after the phase transition the Higgs-field is still ip ~ [EDI. Then, 



since <p < O > h follows from fl2"9|) that m 2 H < 0, and quantum fluctuations grow until the 
evolution of ip becomes potential dominated and the field starts falling to its new minima. 
At this stage the fluctuations of the Higgs-field become impotent and the field configuration 
topology freezes. The fall of ip is very rapid |^0j. After, ip = fx e ff/\/3, the mass of the 
field becomes positive again (|29|). When the field reaches its minimum ip = fi e ff, then, 
m\ = 2A// e// . 
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The magnetic field is formed at the freezing of the Higgs-field configuration. This 
H 2 , |18|], [pL9[| -P| The correlation length £ F at the time of freezing r F , 
or by (153), depending on whether we are still in the linear regime 



\m 



h I 



occurs when 
is either given by 
of not. 



4.2 Computation of £o and £p 

We assume that 4>h is very close to </> , or, equivalently, that the time t h , when the growth 
of £ reaches the speed of light, is very close to the time r = of the phase transition. 
Therefore, 



A = 



< 1 



(35) 



This will be verified when we introduce specific values for the parameters. 
Using (gTI) and <M) we find, 



dr mu 



rripi m 1 



(pa 



Using, (|3~2[) we obtain, 
From (p7[) we have, 
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v/A 7 [<f> 2 H - ^? /2 
2 _ ??vm0o 2/3 



V127TA' 

m pl m 



Thus, solving for th we obtain, 

r H = -v / 3^(127rA , m p ,m0 o ) _1/3 

and thus, 

1 m pl m<pQ 1/3 w m-1/3 

£g = -7= [ /ttt^tt J = 2V37r(127rAm pi m0o) 7 
V A v l^vrA 

Using fl3"3] ) and the above we obtain, 

£o = 3v / 37r(127rA / m p ; m 0o) _1 ^ 3 



(36) 

(37) 
(38) 

(39) 
(40) 

(41) 



Let us compute, now, the correlation length at the time when the field configuration 
freezes, i.e. when \m 2 H \ ~ H 2 . From (|19|), (|23|) and (p9|) we obtain, 



47T m 2 (/)n 



3 A'm^ 



(42) 



2 Nagasawa and Yokoyama |^T| suggest that the freezing of the field occurs a bit later. However, with 
the set of parameters used (see 8.1), the corresponding difference in the correlation length is less than an 
order of magnitude. 
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where 0^ is the magnitude of the inflaton at the time of freezing Tp, for which we find, 



k2 A2 



r F ^^ ^-\ F> „J l ^l^L (43) 
m p im<pQ V 3 \'m 6 pl 

where we have used (27) and the assumption fl35l). 

Using (^), ( |2"3"D and (^) we find that the number of e- foldings of inflation between the 
phase transition and the freezing of the field configuration is given by, 



Itt 2 m 4>q 
"3"AV5 



(44) 



Only during this time are the quantum fluctuations of the Higgs-field important .0 If 
ANp < 1 then the transition proceeds rapidly ||19|| . 

If we are still in the regime of linear growth of the correlation length the value of it at 
the time of freezing is simply given by (|34"D, 



= £o + r F (45) 
If, however, the linear growth of £ has ended before 7>, then it is given by (|28D , 

( 2 ) / 3 rripi 
V o7T mcpo 

where we have used that, m 2 H ~ A'(0q — (/>f) an( ^ (^)- 

Therefore, the initial correlation length of the Higgs-field configuration is given by, 

^ = min(^ 1) ,d 2) ) (47) 



5 Evolution of the magnetic field 
5.1 During inflation 

From (|47|) and ( |14"D we estimate the initial value of the magnetic field produced just after 



the phase transition to be, 

B ~ £" 2 (48) 

Of course, after the GUT phase transition we still have electroweak unification. There- 
fore (H3) represents, in fact, an "electroweak" magnetic field. However, since the residual, 



3 During this time the quantum fluctuations of the Higgs-field could create a "mountain" on the spec- 
trum of the adiabatic density perturbations from the inflaton field (as well as isocurvature fluctuations 
from the Higgs-field) |l8) , [go) . The scale of the perturbations is determined by the Horizon size at the 
time of the transition. If the transition occurs when the observable scales leave the horizon (about 60-45 
e- foldings before the end of inflation), there is a chance of producing excessive density fluctuations and 
microwave background anisotropies. For a detailed review see Salopek et al. [E2|. 
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Figure 1: Evolution of the correlation length £ of the Higgs-field configuration near the 
phase transition. The linear growth of £ starts at r#. £o is the correlation length at the 
time of the transition, ^ is the correlation length at the time the Higgs-field configuration 
freezes, when this occurs during linear growth, and £p is the correlation length at freezing 
when this occurs after the end of the linear regime. 
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electromagnetic symmetry generator is just a projection of the electroweak generators 
(through the Weinberg angle), the residual "electromagnetic" magnetic field will be of the 
same order of magnitude as the one given by (fI8D . Therefore, from now on, we will ignore 
the electroweak transition and treat the above magnetic field as "electromagnetic" . 

During inflation the magnetic field is not frozen into the supercooled plasma |7j but still 
it scales as a~ 2 , since it remains coupled to the Higgs-field ip and, thus, (fH|) is still valid. 
The configuration of the Higgs-field remains comovingly frozen during inflation because 
the initially correlated volumes expand exponentially, faster than causal correlations. This 
is not the case after inflation ends. The scale factor, then, grows slower than the causal 
correlations and the comoving picture of the Higgs-field configuration starts changing as 
the field becomes correlated over larger and larger comoving volumes. However, after the 
end of inflation the magnetic field gets frozen into the reheated plasmafj and decouples 
from the Higgs-field. 

5.2 The rms magnetic field 

In order to estimate the magnetic field on scales larger than the typical dimensions of the 
correlated volumes, we have to introduce a statistical method to do so. 

A thorough treatment by Enqvist and Olesen |12|] suggests that, in all cases that the 
Higgs-field gradients are a diminishing function of n (number of correlated domains), the 
root mean square value of the field would behave as, 

B rms = J(IP) = -^B cd (49) 

where B c d is the field inside a correlated domain and n is the number of correlation length- 
scales, over which the field is averaged. In their treatment Enqvist and Olesen choose the 
Higgs-field gradients as the stochastic variables and also assume that their distribution 
is gaussian and isotropic. Choosing the magnetic field itself as the stochastic variable, 
Enqvist and Olesen reached the same result (^9[) . 

At this point it should be mentioned that in the above treatment the rms value of the 
field has been computed as a line average, that is an average over all the possible curves 
in space between the points that fix the lengthscale, over which the field is averaged. The 
above result may be sensitive to the averaging procedure. One argument in favour of line- 
averaging is that the current galactic magnetic field has been measured using the Faraday 
rotation of light spectra, which is also a line (line of sight) computation. If we assume 
that the ratio of the seed field for the galactic dynamo and the currently observed galactic 
field is independent of the averaging procedure then this would suggest that line-averaging 
is required for the computation of the primordial field. However, the nonlinearity of the 
dynamo process as well as the rather poor knowledge we have for galaxy formation make 
such an assumption non-trivial. In any case, apart from the above, there seem to be 

4 see also section 7.5 
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no other argument in favour of a particular averaging procedure. Therefore, using line- 
averaging could be the safest choice. Here it is important to point out that line averaging 
just gives an estimate of the rms field and does not correspond to any physical process. 

Suppose that we are interested in calculating the rms field at a time t over a physical 
lengthscale L = L(t). Then, 

n(t) = | (50) 

where £ is the correlation length. This scale is equal to the correlation length £o at the 
time of the phase transition. 

In this paper we are mainly interested in the value of the magnetic field at t eq , the 
time of equal matter and radiation densities, when structure formation begins. The scale 
of interest is the typical intergalactic distance, since t eq is preceeding the gravitational 
collapse of the galaxies (see also 7.2). At t eq the corresponding scale is found to be, 

L eq ~ (?f) VS L p ~ 10 pc (51) 

bp 

where t p ~ 10 18 sec is the present time and L p ~ 1 Mpc is taken as the typical intergalactic 
scale at present. 



From ( |50"D the number of correlated domains at t eq is, 

n = n eq = ~p- (52) 

where £ e(? is the correlation length at t eq . 

Therefore, the rms value of the magnetic field over the scale of a protogalaxy is: 

BZs ~ \BZ (53) 



where B e c ^ is the value of the field inside a correlated domain at t eq . 



5.3 Growth of the correlated domains 

It is clear that, in order to calculate the rms field over the galactic scale at t eq we need 
to estimate the correlation length £ e(? , i.e. the size of the correlated domains at that time. 
Therefore, we have to follow carefully the evolution of the correlated domains throughout 
the whole radiation era.Q 

During inflation, as explained already, £ scales as the scale factor a. However, after 
the end of inflation, it grows faster. This is because, when two initially uncorrelated 
neighbouring domains come into causal contact, the magnetic field around the interface is 
expected to untangle and smooth, in order to avoid the creation of energetically unfavoured 

5 The correlated domains should not be pictured as attached bubbles of coherent magnetic field, but as 
regions arround any given point in space in which the orientation of the field is influenced by its orientation 
at this point. 
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magnetic domain walls. In time the field inside both domains "aligns" itself and becomes 
coherent over the total volume. The velocity v, with which such a reorientation occurs, 
is determined by the plasma, which carries the field and has to reorientate its motion for 
that purpose.f] 

Thus, the evolution of the correlation length is given by, 

^ = H£ + v (54) 

where v is the peculiar, bulk velocity, determined, in principle, by the state of the plasma. 

From fl54|) it is apparent that the correlated domains could grow faster than the Hubble 
expansion. Therefore, the magnetic field configuration is not necessarily comovingly frozen 
and the domains could expand much faster than the universe, resulting in large correlations 
of the field and high coherency. 

In order to describe the evolution of the correlated domains one has to determine the 
peculiar velocity v. This primarily depends on the opacity of the plasma. 

If the plasma is opaque on the scale of a correlated domain, then radiation cannot 
penetrate this scale and is blocked inside the plasma volume. Consequently, the plasma is 
subject to the total magnetic pressure of the magnetic field gradient energy. Therefore, this 
energy dissipates through coherent magnetohydrodynamic oscillations, i.e. Alfven waves. 
In this case, the peculiar velocity of the magnetic field reorientation, is the well known 
Alfven velocity 

v A = ^ (55) 

where B c ^ is the magnitude of the magnetic field inside a correlated domain and p is the 
total energy density of the universe, since, before t eq , matter and radiation are strongly 
coupled.0 

If the plasma is not opaque over the scale £ of a correlated domain, then radiation can 
penetrate this scale and carry away momentum, extracted from the plasma through Thom- 
son scattering of the photons. This subtraction of momentum is equivalent to an effective 
drag force, F ~ po~TVT n e @- Balancing this force with the magnetic force determines the 
"Thomson" velocity over the scale £, 

v 2 

v T = (56) 

where va is the Alfven velocity, n e is the electron number density and o~t is the Thomson 
cross-section. 



6 Note that the plasma does not have to be carried from one domain to another or get somehow mixed. 
Also, conservation of flux is not violated with the field's rearrangements, since the field always remains 
frozen into the plasma, which is carried along. 

7 Unless explicitly specified, natural units are being used (h = c = 1). In natural units, G — m~ ; . 

8 This coupling implies that any reorientation of the momentum of matter has to drag radiation along 
with it. This increases the inertia of the plasma, that balances the magnetic pressure. 
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Hence, for a non-opaque plasma the peculiar velocity of the plasma reorientation is 
given by @, 



v = min(t>A, vt) (57) 

In order to explore the behaviour of the opaqueness of the plasma, we need to com- 
pare the mean free path of the photon l T ~ (n e crr) _1 to the scale £ of the correlated 
domains. For realistic models, the correlated domains remain opaque at least until the 
epoch t an h ~ 0.1 sec of electron-positron annihilation (T ~ lMeV). The reason for this 
can be easily understood by calculating It before and after pair annihilation. 

For T > 1 MeV, instead of the usual Thomson cross-section or, we have the Klein- 
Nishina cross-section 



3 ,m Ps r , IT 1 . .GeV , , r T , ^ , . 

- g ( j?) [In - + - ] = 2.7 (_) !»[_] GeV (58) 

where m e ~ 0.5GelTs the electron mass and ot — 6.65 x 10~ 25 cm 2 ~ 1707.8 Gey -2 . The 
electron number density is given by p5| , 



n e ^- A C -^ 9e T* (59) 

where £(3) ~ 1.20206 and g e = 4 are the internal degrees of freedom of electrons and 
positrons. 

From (|58~D and (|59|) we find, 

It ~ 0,1 il er for T > 1 MeV (60) 

T l 

which at annihilation gives, lT{t an h) ~ lO^GeV -1 . 

After annihilation the electron number density is given by ||25|| , 

n e ~ 6 x 10~ 10 n 7 ~ 1.44 x 10~ 10 T 3 (61) 
where n 7 is the photon number density given by, 

n ^M g T 3 (62) 

where g 1 = 2 are internal degrees of freedom of the photon. 
With the usual value for <jt we obtain, 

10 6 GeV 2 

It ~ ^ for T < 1 MeV (63) 

At annihilation the above gives, lT(t a nh) ~ 10 15 GeV . 

Hence, the mean free path of the photon at the time of pair annihilation is enlarged by 
a factor of 10 10 ! As a result, It is very likely to become larger than £ after t an h- If this is 
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so, the Thomson dragging effect has to be taken into account and the peculiar velocity of 
the plasma reorientation is given by (|57|) . 

In order to calculate the peculiar velocity it is necessary to compute the Alfven velocity, 
which requires the knowledge of the magnetic field value B cd inside a correlated domain. 
To estimate this we assume that the magnetic flux, on scales larger than the sizes of the 
correlated domains, is conserved, as implied by the frozen-in condition. 

Consider a closed curve C in space, of lengthscale L > £, encircling an area A. Conser- 
vation of flux suggests that the flux averaged mean magnetic field inside A scales as a~ 2 . 
This implies that for the field inside a correlated domain we have, B c d(L/ 'C,)^ 1 oc a~ 2 .f\ 
Since C follows the universe expansion L oc a, with a oc t 1 / 2 . Thus, for the radiation era, 
we obtain, 

B cd t l ' 2 i = K^B cd =^ (64) 

where K is a constant to be evaluated at any convenient time. Since, the correlation length 
grows at least as fast as the universe expands, the magnetic field inside a correlated domain 
dilutes at least as rapidly as a~ 2 for the radiation era. 
Substituting the above into (BH) we find, 



K t 1/2 

va-IQ (65) 

mpi 4 

Solving the evolution equation (|54] ) with a oc t 1 / 2 in the case that v — Va gives, 

m 2 = {j)& +^A(*K(t)t(i-yf ) (66) 

where & is the correlation length of the field at the time ti. The first term of (|66|) is due 
to the Hubble expansion, whereas the second term is due to the peculiar velocity. 

In the case of v — vt, for t > t an h, using fl6~i"]), and the usual value of <tt, fl56|) gives, 



f 5/2 

where 



v T = D — (67) 



D ~ 10- 57 K 2 GeV~ 3/2 (68) 
Using flBTD, the evolution equation ( p^ ) gives, 

^) 4 = (^) 2 C 4 + ^TW^[l-(f) 3/2 ] (69) 

9 Note that the flux averaging of the field on scales larger than the correlation length, corresponds to 
a physical process, that of the field untangling, and is so in order to preserve flux conservation on scales 
that the magnetic field is frozen into the plasma. This should not be confused with the line averaging 
procedure which we use to estimate the rms field, and does not correspond to a physical process. 



15 



The evolution of the correlation length of the magnetic field configuration is described 
initially by the Alfven expansion equation ( |S6f) until the moment when £ ~ Zy. From then 
on the growth of £ continues according either to (|66|) or to fl69|), depending on the relative 
magnitudes of the velocities va and vt- Using the above, we can calculate the scale £ eg of 
the correlated domains at t eq and, thus, calculate the rms magnetic field from (|53|). 



5.4 Diffusion 

An important issue, which should be considered, is the diffusion length of the freezing of 
the field. Indeed, the assumption that the field is frozen into the plasma corresponds to 
neglecting the diffusion term of the magnetohydrodynamical induction equation [p6|j , 



— = V x (v x B) + <r -1 V 2 B (70) 
at 

where v is the plasma velocity and a is the conductivity. In the limit of infinite conductivity 
the diffusion term of ( |T0D vanishes and the field is frozen into the plasma on all scales. 
However, if a is finite then spatial variations of the magnetic field of lengthscale I will 
decay in a diffusion time, r ~ al 2 Thus, the field at a given time t can be considered 
frozen into the plasma only over the diffusion scale, 

U ~ (71) 

If Id > £, the magnetic field configuration is expected, in less than a Hubble time, 
to become smooth on scales smaller than ld(t). Thus, in this case, it is more realistic 
to consider a field configuration with coherence length and magnitude of the coherent 
magnetic field B c d, where B c d = B^/n^ is the flux-averaged initial magnetic field over 
n d = ld/ii number of domains. 

An estimate of the plasma conductivity is necessary to determine the diffusion length. 
The current density in the plasma is given by, J = nev, where n is the number density 
of the charged particles. The velocity v acquired by the particles due to the electic field 
E, can be estimated as v ~ eEr c /m, where m is the particle mass and r c = l m jp/v is the 
timescale of collisions. Since the mean free path of the particles is given by, l m fp ~ l/na c , 
the current density is, J ~ e 2 E/mva c , where a c is the collision cross-section of the plasma 
particles. Comparing with Ohm's law gives for the conductivity ||26|| , [27], 



e 2 



a 



mvop 



(72) 



The collision cross-section is given by the Coulomb formula |27_ 



^e^lnA (73) 

where In A ~ \n(e~ 3 ^jT 3 /n) is the Coulomb logarithm. Thus, the behaviour of the conduc- 
tivity depends crucially on the temperature. 
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For low temperatures, T < m e ~ 1 MeV (i.e. after t an h), the velocity of the electrons 



is, v ~ yT/m e . Thus, from (|T2| ) and (|7^) the conductivity is given by, 

(74) 




m e In A 

For high temperatures, T 3> m e , (^J) suggests that, In A ~ 1. Also, the mass of the 
plasma particles is dominated by thermal corrections, i.e. m ~ T, and v ~ 1. Consequently, 
in this case, ( |72|) and ([T3|) give for the conductivity, 

o ~ ^ (75) 

Using the above results we can estimate the diffusion length. Indeed, from ([711) , ( |74| ) 
and ( |75| ) we obtain, 

f 10 8 GeV 1/2 T- 3 / 2 T > 1 MeK 

Id ~ ^ (76) 
[ 10 8 Ge1/ 3 / 4 T- 7 / 4 T < 1 MeV 

An important point to stress is that the diffusion length is also increasing with time. If 
Id > £, then the size of the correlated domains is actually determined by the diffusion length 
and it is the growth of the later that drives the evolution of the magnetic field configuration. 

At this point we could briefly discuss the behaviour of any electric field, produced by the 
phase transition. As can be seen by (|5|) and ([T2|), the electric field, = is determined 
by the time derivative of the Higgs-field VEV. Thus, strong currents are expected to arise 
at the time of the transition, when the VEV of the Higgs-field falls rapidly from zero to fi. 
These are the currents that accompany the creation of the magnetic field [f29fl . However, 
after the transition and during inflation, the VEV of the Higgs-field, at any point in space, 
is more or less fixed and constant in time, since the field configuration is comovingly frozen. 
Therefore, there should not be any significant electric field surviving the transition. After 
inflation this comoving picture begins to change but the magnetic field decouples from the 
Higgs-field and, thus, any electric field produced, by shifting of the magnetic field lines, 
is related to plasma motion phenomena. Since such reorientations occur, we expect small 
electric fields to be present in the form of electromagnetic waves, which will diffuse and 
thermalize the gradient energy of the magnetic field, that is reduced by its reorientation 
and alignment. 

6 At the end of inflation 

6.1 The Reheating Temperature 

The time t en d when inflation ends could be determined by the reheating temperature T reh 
with the use of the well known relation, 

W^0.3s; 1/2 (^) (77) 

reh 
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where g* is the number of particle degrees of freedom, which, in most models, is of order 10 2 
(e.g. in the standard model it is 106.75 whereas in the minimal super symmetric standard 
model is 229). 

The reheating temperature is usually estimated by f30fl, 



^ « 0.78 a VV 1/4 (— ) 1/a (7g) 
m p i m p i 

where H en d is the Hubble parameter at the end of inflation and a is the reheating efficiency, 
which determines how much of the inflaton's energy is going to be thermalised. Using ( |T9"D 
and ([24]) we obtain, 

Hend ^ 77= (79) 



Subsituting to (|77|) and (|78"D we find, 

T 2 reh ^ 0.35 a^g^mmpt (80) 

In most inflationary models reheating is prompt, it is completed quickly and a ~ 1. In 
case of a quadratic inflaton potential, however, as in false vacuum inflation, the reheating 



process could be incomplete and extremely inefficient |3T[. However, the magnitude of 
the reheating inefficiency is still an open question. Kofman et al. [3^] suggest that the 
reheating temperature would be of the order, T reh ~ lO^^/mwy which, compared to (|80|), 
implies that a ~ 10~ 4 . Shtanov, Trashen and Brandenberger |H| make a lower estimate, 
Treh ~ m. 

6.2 Thermal fluctuations 

The Higgs-field, through the Higgs mechanism, provides the masses of the particles after 
the GUT phase transition. Thus, it is in that way coupled to the thermal bath of the 
particles. Therefore, at reheating, this coupling introduces thermal corrections to the 
effective potential of the Higgs-field. Consequently, if the reheating temperature is high 
enough, the configuration of the Higgs-field may be destroyed due to excessive thermal 
fluctuations. This will eraze any magnetic field if the later has not been frozen into the 
plasma already. The above will occur if the temperature exceeds the well known Ginzburg 
temperature Tq- Moreover, if the temperature exceeds a critical value T c there is danger 
of thermal restoration of the GUT symmetry itself. 



The Ginzburg and the critical temperatures are simply related |34 

T C -T G ~\T C (81) 

Thus, for A < 1, 

T G ~T C ~ y/Xfx (82) 

Therefore, it is very important to see if the temperatures during the reheating process 
could exceed T c . At this point it should be noted that the reheating temperature is not 
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the highest temperature achieved during the reheating process. Indeed, as soon as the field 
begins its coherent oscillations, the temperature rises rapidly and assumes its maximum 



value M, 25 



Tmax — (yendTreh) 1 ^ (83) 



where V en d is the energy density of the inflaton at the end of inflation. From ( p3|) and 
we obtain, 

V end ~ 0.1 m 2 m 2 pl (84) 

Thus, 

T m ax ~ {yjmm p i T reh ) 1/2 (85) 

Therefore, in order to avoid symmetry restoration and any Ginzburg phenomena we should 
have, 

X/i 2 

Tmax < T c =>■ T ren < (86) 

If the reheating temperature exceeds the above value then the magnetic field is ther- 
mally unstable and we are in danger of restoring the GUT symmetry. However, if the field 



survives then its stability is ensured |29 



After reaching its highest value, T max , the temperature slowly decreases during the 
matter dominated era of the coherent inflaton oscillations, until it falls to the value T ren 
when the universe becomes radiation dominated.PI 



7 Constraints 

7.1 Constraints on the parameters 

If we assume that the observed density perturbations are due to inflation, then we have 
from COBE 

v N7r m = 5.5 x 1(T 6 (87) 



m p i 

which yields, 

m ~ 10 l3 GeV 

Other restrictions of the model imposed on /i, A and A' are |2e 

0<A,A'<1, /I < ^= ~ 10 18 GeV (89) 



57T 

An additional constraints for the A are established by the inflaton domination condition 
j|). Also the ratio of the A's can be determined by ([H]) and ( p6|) with the reasonable 



10 We should mention that this small period of matter domination is not taken into account in our 
treatment due to the fact that its duration is very small compared to the timescales considered and so we 
choose to ignore it for the sake of simplicity. 
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assumption that the number of e-foldings of inflation after the phase transition is of order 
N ~ 10, 

p ~ (90) 
A' p 

7.2 The Galactic Dynamo constraint 

From the present understanding of the galactic dynamo process || it follows that, in order 
for a primordial magnetic field to be the seed for the currently observed galactic magnetic 
field, it should be stronger than 10~ 19 Gauss at the time of galaxy formation, on a comoving 
scale of a protogalaxy (~ 100 kpc). 

Since the gravitational collapse of the protogalaxies enchances their frozen-in magnetic 
field by a factor of (pc/Pc) 2 ^ 3 ~ 10 3 (where pc ~ I0~ 24 gcm~ 3 is the typical mass density 
of a galaxy and p c ~2x lO~ 29 Qh 2 g cm~ 3 is the current cosmic mass density), the above 
seed field corresponds to an field of the order of ~ 10~ 22 Gauss over the comoving scale of 
~ 1 Mpc. With the assumption that the rms field scales as a~ 2 with the expansion of the 
universe (a oc t 2 ^ 3 for the matter era), we find that the required magnitude of the seed field 
at t eq is ~ lO~ 22 Gauss x (t gc /t eq )^ 3 ~ lO~ 20 Gauss, where t gc ~ 10 15 sec is the time of the 
gravitational collapse of the galaxies. 

The above justify our choice to calculate the magnetic field at t eq over the comoving 
scale of 1 Mpc and consider the constraint, 



B eq > 10~ 20 Gauss (91) 

From recombination onwards the non-linear nature of structure formation is very dif- 
ficult to follow. Indeed, there exist a numerous collection of different models. A strong 
primordial magnetic field could influence in various ways some of these models, possibly 
with a positive rather than a negative effect. 



7.3 The Nucleosynthesis constraint 

One upper bound to be placed on the magnetic field at t eq is coming from nucleosynthesis. 
This has been studied in detail by Cheng et al. |55|]. They conclude that, at t nuc ~ 1 sec, 
the magnetic field should not be stronger than, 



B nuc < 1Q n Gauss 



(92) 
relaxes 



on a scale larger than ~ 10 4 cm. A more recent treatment by Kernan et al. [^G 
the bound by about an order of magnitude, B nuc < e~ x {T™ uc ) 2 ~ 10 12 Gauss, where T u is 
the neutrino temperature and e is the electric charge. This bound is valid over all scales. 
Similar results are also reached by Grasso and Rubinstein J37 ] . 
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7.4 Energy Density constraints 

Constraints are also induced by ensuring that the energy density of the magnetic field is 
less than the energy density of the universe. During the inflationary period, due to inflaton 
domination, the energy density of the universe is mainly in the inflaton field. However, 
after reheating and until t eq , the energy density of the Universe is just the radiation energy 
density. 

Thus, for the inflationary period we should verify that, 

l>n <1 (93) 



Pinf 

where ps = B 2 d /8n and p in j = V(<p) are the energy densities of the magnetic and inflaton 
fields respectively. 

The highest value of the above ratio corresponds to the time of the phase transition 
since the magnetic field is rapidly diluted during inflation, whereas the inflaton's potential 
energy remains almost unchanged. Using fl4*8] ) and (^) we find the first energy density 
constraint, 



m0 o > & (94) 

After reheating, the expansion of the universe dilutes the energy density pb of the 
magnetic field, inside a correlated domain more effectively than the radiation density, 
which scales as a -4 . Therefore, it is sufficient to ensure that ps(t) is less than the energy 
density p(t) of radiation at the time t« of the formation of the magnetic field configuration. 
That is, 

< =*3* < ^ ^ (95) 
which is the second energy density constraint. 

7.5 The non-Abelian constraint 

During the electroweak era, the freezing of the magnetic field into the electroweak plasma is 
not at all trivial to assume. Indeed, before the electroweak transition, since the electroweak 
symmetry group SU(2) x U(l)y is still unbroken, there are four apparent "magnetic" fields, 
three of which are non-Abelian. 

It would be more precise, then, to refer only to the Abelian (Hypercharge) part of the 
magnetic field, which satisfies the same magnetohydrodynamical equations as the Maxwell 
field of electromagnetism. The non-Abelian part of the field may not influence the motion of 
the plasma due to the existence of a temperature dependent magnetic mass, m B « 0.28 g 2 T 



(see for example |38| and P9"|), which could screen the field over the relevant lengthscales. 

The condition for this screening to be effective can be obtained by comparing the 
screening length r$ ~ m^ 1 of the non-Abelian magnetic fields with the Larmor radius of 
the plasma motion rx ~ jg, where m ~ \faT (a = g 2 /4ir) is the temperature induced 
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physical mass of the plasma particles, g ~ 0.3 is the gauge coupling (charge) and v is the 
plasma particle velocity. If we assume thermal velocity distribution, i.e. mv 2 ~ T, we find, 

r , t 2 

R = — ~ 10~ 2 — (96) 
rs B cd 

If R > 1 then our restriction to the Abelian (Hypercharge) part of the magnetic field 
is well justified. This restriction will not cause any significant change to our results since, 
at the electroweak transition, the hypercharge field projects onto the photon through the 
Weinberg angle (|), cos 6w ~ 0.88. If, however, R < 1 then the non- Abelian fields do affect 
the plasma motion, and should be taken into account. Since T oc a^ 1 and B c d falls at least 
as rapid as a~ 2 , R is, in general, an increasing function of time. Thus, the constraint has 
to be evaluated at reheating. 



7.6 The monopole constraints 

Unfortunately, the mechanism, which we use to generate the primordial magnetic field, 
could also produce stable magnetic monopoles. Since these monopoles should not dominate 
the energy density of the universe, we require that the the fraction fl M of the critical density, 



contributed by the monopoles, to be less than unity, that is E5 



n-^^THlo^* 1 (97) 

where M = Airfig' 1 ~ 10 [i is the monopole mass, um is the monopole number density, s is 
the entropy density of the universe and h is the Hubble constant in units of 100 km/sec/Mpc. 
The ratio Um / s is a constant^ and can be evaluated at the end of inflation t en d- Taking 
n M ~ iendi wnere iend is the correlation length at that time, we have, 

UM^ TreH, 
S C m pl 

where ( = £, e ndH en d gives the correlation length as a fraction of the Hubble radius. From 
(p7|) and ( p8|) we find the first monopole constraint p3| , 



C 3 >10 n ( Treh f( - ) (99) 

^ ~ l 10 14 GeV rJ ho^GeV' 1 ' 

Apart from the above, mass density constraint, another constraint is the well known 



"Parker bound" [40], which considers the effect of galactic magnetic fields onto the magnetic 



monopole motion. The flux of the monopoles is 



$m = ^n M v M ~ lO^X-^cm-V-W 1 (100) 

47T S 10 -ci 



1 we can ignore monopole annihilations (see |25| ) 
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where % is the monopole velocity. The monopoles are accelerated by the galactic magnetic 
field B g ~ 10~ 6 Gauss to velocity, 

2huBl l/2 ^ 3 10^GeV 1/2 

M M y ! 

where I ~ 1 kpc is the coherence length of the magnetic field and hu ~ e _1 is the magnetic 
charge of the monopole. 

The magnetic field ejects the monopoles from the galaxy, while providing them with 
kinetic energy Ek ~ hB g l ~ 10 u GeV. Demanding that the monopoles do not drain the 
field energy in shorter times than the dynamo timescale, i.e. the galactic rotation period 
t ~ 10 8 yrs, we find the constraint, 

B 2 II 

g/ > t => $ M < 10~ 15 cm- 2 sr' 1 sec' 1 (102) 



where ~ 30 kpc is the size of the galactic magnetic field region. Using and from 
( |100| ) and ( |102j ) we find the second monopole constraint, 

c3£1 ° 12( io^F )3(12 T^ )1/2 (103) 

7.7 Additional constraints and considerations 

Finally, we have to make sure that at the time the magnetic field is formed the correlation 
length given by (|47D is still inside the horizon, that is, 

Hp 1 > & (104) 

where Hp is the Hubble parameter at the time of the formation of the magnetic field. 
With great accuracy Hp ~ Hq, where Hq is the Hubble parameter at the time of the phase 
transition. H can be easily computed using (p~6f) , (|20"D, and (p7|) for the potential (|23|). 



8 Evaluating 

In order to be consistent with our assumptions (e.g. inflaton domination) we will consider 
the phase transition to take place at the latest at, N ~ 1. We first choose a set of typical 
model-parameters. 

8.1 Choosing the values of the parameters 

As already mentioned, the mass of the inflaton field m is determined by COBE, 

m ~ 10 13 GW (105) 
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For the selfcoupling of the Higgs-field A we choose the usual value, 

A ~ 1 (106) 

Inserting the above values into the inflaton domination condition ( |25| ) we find that the 
maximum value of \x is, 

fi ~ 10 15 GeV (107) 

Finally, the coupling A' between the Higgs-field and the inflaton can be determined 
with the use of Q90|) , 

A' ~ 10~ 8 (108) 



8.2 For all N 



Now that the parameters of the model are chosen, the only parameter still to be determined 
is the number N of e-foldings of inflation, which remain after the phase transition. We will 
treat this as a free parameter, link it with the resulting magnetic field and then try, with 
the use of the constraints, to establish its extreme values. In that way we will be able to 
fully examine the corresponding behaviour of the field at t eq . 

We begin by extracting some direct, iV-independent results from the, previously chosen, 
values of the model-parameters. From ( |S6"D we find that the upper bound for the reheating 
temperature is estimated to be, 

T reh ~ 10 14 Ge1/ (109) 

which is in agreement with the estimates of Kofman et al ]32| and higher than the estimates 
of Shtanov et al [33]. 

Using flTTD this gives the time when inflation ends, 



With the use of flTBD and 



tend ~ 10- 35 SeC 

the assumption (^f 
A ~ 10~ 2 < 1 



D is easily verified, 



;no) 



(in; 



Now, for the correlation length, from plf) we find, 



(112) 



From (16), and 



Thus, from 



we obtain, 



Using ( 46|) we also find, 



10- Lb GeV- L < £ 



(i) 



W~ u GeV- 



.(2) 



10~ u GeV- 



(113) 

(114) 
(115) 
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Therefore, from ( |i?D and the above we have, 

£ F ~ £ ~ lO^GeV- 1 ~ 10" 28 cm ~ 10" 4 V (116) 
Now, from (fHj) we find that, 

AiV F ~ 10~ 3 < 1 (117) 

and, therefore, the phase transition is very rapid. 

We can, now, check on the horizon constraint (|104 ). The value of Hp is found to be, 

H F ~ Hq ~ 10 13 GeV (118) 

Comparing with (|116|) we see that the constraint is satisfied. By using ( |116|) into ( 94]) we 
find, that the first energy density constraint is also satisfied. 
The initial magnetic field is found from ( |48|) and (|116|) to be, 

B ~ 10 47 Gauss (119) 

Let us now evaluate the A^-dependent quantities. 
The correlation length at the end of inflation is, 

Ud = —Zf ~ lO-^GeV- 1 (120) 
a 

where we have used (|22|) , (|26|) and ( |116| ). 

From the above and considering also the fact that, during inflation, the magnetic field 
configuration is comovingly frozen, we find that the magnitude of the magnetic field inside 
a correlated domain is given by, 

B e c f ~ 10 4 V 2A 'Gauss (121) 
Evaluating, (|96|) at the end of inflation, we find, 

i?~0.1e 2Ar >l foriV>l (122) 

and the non-Abelian constraint is satisfied for all N. 

Using ( |110| ) and ( p.21| ) we can show from ( Po]) that the second energy density constraint 
is also satisfied for all N. 

From (|3|) we find, 

K ~ l0 8 e~ N GeV 1/2 (123) 

We evaluated the above also at the end of inflation, using (|H|), ( |1 10| ) and ( |120| ). 

At early times the correlated domains are opaque to radiation and, thus their growth 
is determined by ( pB| ) with tj — > t en <i. The domains remain opaque at least until the time 
of the electron pair annihilation 

At annihilation t an h ~ 0.1 sec, (|B"6|) gives, 
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10 3 e N GeV- 1 N > 15 

Uh ~ <( (124) 
10 i3 e -7v/2 Ge y-i l < < 15 

Comparing with the photon mean free path it is evident that, £ an /j > It ~ 10 5 GeV A_1 , 
for all iV. 

However, after annihilation, It increases drastically in size, lT{t a nh) ~ lO^GeV -1 . Com- 
paring this value with ( |124] ) we find that, at T ~ 1 MeV, the correlated domains become 
transparent to radiation for N < 28. Thus, for N >28 the Alfven expansion continues 
after pair annihilation, whereas, for N < 28 the Thomson scattering effect has to be taken 
into account. 

CASE 1: For 1 < N < 28 

If N < 28 then, after t an h the Thomson effect has to be taken into account. 
The Alfven velocity va at t an h is found by (|65|) with the use of ( |123|) , 

0.le~ 2N A>15 

v A (t anh ) ~ { (125) 
10 -ii e -iv/2 1< AT < 15 



Similarly, the Thomson velocity vt at t an h is found by (^3) with the use of fl64|) and 

f 10 10 e- 5Ar > 15 

MWO ~ i (126) 
{ lO- 20 e' N / 2 1 < N < 15 

From the above it is straightforward that, for N < 28, VT(t an h) < VA(t an h). Therefore, 
after t an h, the evolution of the correlated domains is determined by the Thomson effect. If 
we assume that the Alfven expansion does not take over again until t eq then the correlation 
length at that time can be obtained by (|69|), 



10 9 e N GeV- 1 19 < N < 28 

U ~ { (127) 
10 2i e -A72 Ge y-i 1 < AT < 19 

Using this value we can verify that the Thomson velocity remains always smaller than 
the Alfven velocity until t eq . Physically, (|127|) implies that, if N > 19 the damping of the 
growth of the correlated domains, after t an h is so effective that the Hubble term dominates 
their evolution. However, for 1 < N < 19, the Thomson velocity, although small, is still 
capable of outshining the Hubble term. 

CASE 2: For N > 28 

For high values of N the magnetic field is so much diluted by inflation that the Alfven 
or Thomson expansions are insignificant. The growth of the correlated domains is driven 
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solely by the Hubble expansion and, thus, 



i eq ^\h^~ Uh ~ !0 9 e N GeV- 1 N > 28 (128) 

V ^anh 

In total, ( p.27| ) and Ql28|) suggest the following behaviour for the correlation length at 

( 10 9 e N GeV- 1 N > 19 

U ~ (129) 
I KFe-^GW- 1 1 < N < 19 

From (|76| ) we find that at t eq the diffusion length is, l^ 9 ~ 10 23 £reV -1 . Comparing with 
the above we see that l e £ > £ e <j for N < 32. Thus, the dimensions of the correlated domains 
at t eq are actually given by, 

f 10 9 e N GeV~ 1 N > 32 

U ~ (130) 
[ ^ 9 ~ W^GeV' 1 1 < iV < 32 

8.3 The magnetic field's range of values 

We are, now in the position to calculate the magnetic field strength at t eq . From fl64]) and 



10 e- m Gauss N > 32 

Kl ~ { (131) 
10- 13 e~ N Gauss 1 < iV < 32 

Also, from ([52] ) and (|130j) we get, 

f l0 24 e~ N N > 32 

7i ~ < (132) 
[ 10 10 1 < N < 32 

With the use of the above, in view also of (fxf), we can immediately find the rms value 
of the field at t eq for a given N, 

lO-n e -3N/2 Gauss N > 32 

BZs ~ <! (133) 
10 -i8 e -JV GaMSS 1< iV < 32 



|130| ) we have, 



As can be seen from (|133|) , the maximum rms value of the field at t eq corresponds to 
N = 1, 

(BZs)ma X ~ lO- ls Gauss (134) 
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For the minimum value of the field we just employ the galactic dynamo constraint (|9l|). 
This gives, 

N max ~ 5 (135) 
Thus, the range of values of the magnetic field is, 

5 > N > 1 
10~ 20 Gauss < B £ i < 10- 18 Gauss (136) 

The above result, however, are still subject to the nucleosynthesis and monopole con- 
straints. 

8.4 Nucleosynthesis and monopole constraints on N 

Since nucleosynthesis occurs very near the electron pair annihilation we will assume, for 
simplicity, that the correlation length £ nuc at t nuc ~ 1 sec is approximately equal to the one 
at annihilation, i.e. £ nttc ~ £ anh JQ 

The diffusion length at t anh is found by fl76|) to be, l a d nh ~ lO^GeV -1 . Thus, from ( ggg) 
we have, 

r lO^^^GeV- 1 N > 23 

£nite ~ ianh ~ S (137) 

[ ^ lO^GeV- 1 1 < iV < 23 
Inserting the above into (|63|) and with the use of (|123|) we find, 

10 13 e~ 2N Gauss N > 23 

(138) 

10 3 Gawss 1 < < 23 

The maximum value of the magnetic field is, B'^ C (N — 1) ~ 10 2 Gauss. Comparing 
with (P2]), we see that the maximum value of the field is well below the nucleosynthesis 
constraint and, therefore, the constraint is not violated for any value of N. 

Let us consider the monopole constraints. Given the values of the model parameters 
and the assumed reheating temperature, both of the monopole constraints (|99| ) and ( |103| ) 
reduce to, C > 10 4 . Using (|79D and fll20D we find that, 

C = UdH end ~Q.le N (139) 

and the constraints are satisfied only if AT > 11. Thus, a magnetic field strong enough to 
seed the galactic dynamo, would violate the monopole constraints. 

One way to overcome the monopole problem is to consider GUT-models which do 
not admit monopole solutions, such as "flipped" SU(5), i.e. the semi-simple group, 
SU(5) x £7(1). 

12 The magnitude of £„ uc does not affect the results when the Hubble term in the evolution equations is 
subdominant. In the opposite case our assumption perturbs the results by an order of magnitude in the 
values of £ eq and n but less than an order of magnitude in the value of the rms magnetic field, since the 
later depends on 1/y/n [see (p9[)]. 



B cd ~ 
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9 Conclusions 



We have analysed the creation and evolution of a primordial magnetic field in false vac- 
uum inflation. We have shown that, in GUT-theories that do not produce monopoles, 
a sufficiently strong primordial magnetic field can be generated, provided that the phase 
transition takes place no earlier than 5 e-foldings before the end of inflation. Although the 
magnetic field produced is strong enough to seed the dynamo process in galaxies, it does 
not violate any of the numerous constraints imposed (apart from the monopole constraint 
if applicable). 

Our results are sensitive to the reheating efficiency. Indeed, if reheating is efficient, then 
the time of the end of inflation is earlier and the resulting field diluted by the expansion of 
the universe. More importantly, though, if the reheating temperature is of the order of the 
critical temperature or the Ginzburg temperature, then the magnetic field will be erazed. 
Fortunatelly, this does not appear to be the case. 

Finally, the strength of the magnetic field produced by our mechanism relies on the 
value of N, i.e. on the moment that the phase transition occurs. In turn, this depends on 
the exact values of the model-parameters. Observational data on the primordial magnetic 
field could determine, or constraint, these parameters. Experiments to detect such a field 



have occasionally been suggested (see for example |41| or f42|). Not merely would the 
observation of a primordial field yield information on false vacuum inflation, but it would 
also improve our understanding of the galactic dynamo and of non-linear astrophysical 
processes in general. 
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